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The long time accumulation of the random actions of a single particle "reservoir" on its coupled system can 
transfer some temperature information of its initial state to the coupled system. This dynamic process can be 
referred to as a quantum thermalization in the sense that the coupled system can reach a stable thermal equi- 
librium with a temperature equal to that of the reservoir. We illustrate this idea based on the usual micromaser 
model, in which a series of initially prepared two-level atoms randomly pass through an electromagnetic cavity. 
It is found that, when the randomly injected atoms are initially prepared in a thermal equilibrium state with a 
given temperature, the cavity field will reach a thermal equilibrium state with the same temperature as that of the 
injected atoms. As in two limit cases, the cavity field can be cooled and "coherently heated" as a maser process, 
respectively, when the injected atoms are initially prepared in ground and excited states. Especially, when the 
atoms in equilibrium are driven to possess some coherence, the cavity field may reach a higher temperature in 
comparison with the injected atoms. We also point out a possible experimental test for our theoretical prediction 
based on a superconducting circuit QED system. 

PACS numbers: 05.30.-d, 03.65.Yz, 85.25.-j 



I. INTRODUCTION 

A small system in contact with a large reservoir (or so- 
called heat bath) in thermal equilibrium of temperature T will 
dynamically approach to an equilibrium state with the same 
temperature T QJJ] . This irreversible process from a nonequi- 
librium state into a stable one is conventionally referred to as 
quantum thermalization. Most recently, another kind of ther- 
malization, called canonical thermalization (e.g., Refs. 10-Hll), 
investigated in the meaning of typicality that almost all pure 
states in the universe (the system plus its bath) are entangled, 
and thus the system can reach an approximately canonical 
thermal state by averaging over the bath. Here, the temper- 
ature appears as an "emergent" concept. 

In conventional thermalization, the heat bath consists of a 
very large number of degrees of freedom (for example, a set of 
harmonic oscillators for the bosonic heat bath), and the cou- 
pling strengths of the thermalized system with the degrees of 
freedom of its bath are randomly distributed. According to 
the viewpoint in statistical mechanics that an average over an 
ensemble is equivalent to the time average in some sense [6], 
a natural question is if a series of random actions of a single- 
particle "reservoir" injected randomly in a time domain can 
transfer some temperature information of its initial state to the 
coupled system at a steady state as a thermalization process? 
To answer this question, in this paper we study the steady state 
of a quantum system which is controlled to have a randomly 
"multipulse" type interaction with a single-particle system ini- 
tially prepared in thermal equilibrium with a temperature. If 
the steady state of the quantum state is a thermal one with 
the same temperature as that of the single-particle system, we 
think that this quantum system has been thermalized by the 
single-particle system through a randomly "multipulse" type 
interaction. 

Since the randomly "multipulse" type interaction can be re- 
alized by random injections, in this paper we will illustrate our 
idea based on the usual micromaser model (e.g., Refs. I7l-l20ll). 



in which a series of initially prepared atoms pass through an 
electromagnetic cavity. Here, the single-mode cavity field is 
the system to be thermalized and the randomly injected atoms 
play the role of the single-particle reservoir 02111 . Under some 
conditions we will clarify, if the injected atoms is initially pre- 
pared in thermal equilibrium, that the conventional thermal- 
ization enables the cavity field to transit from any initial state 
to a thermal state with the same temperature as that of the 
atoms. We also find that the temperature of the cavity field in 
thermal equilibrium depends on the initial state of the injected 
atoms. As in two limit cases, such quantum thermalization can 
describe the cooling [22] and masering processes 0-53 E3l> 
which respectively correspond to the cases where the injected 
atoms are initially prepared in ground and excited states. 

It is worth noting that when the atoms initially possess some 
quantum coherence ll24l l25tl . the "thermalized state" of the 
cavity field will carry the information of this coherence. Ac- 
tually, quantum coherence has been proved to be a kind of 
resource to enhance quantum information processing. Most 
recently, some studies have shown that physical processes 
with quantum coherence usually possess some novel effect 
for energy transfer 11261 12711 . For example, quantum heat en- 
gines using quantum matter (and even with the assistance of 
Maxwell's demon) as a working substance can improve work 
extraction as well as the working efficiency in the thermody- 
namics cycle 112 8143 Oil . In the present study, it is expected that, 
when the injected two-level atoms possess some coherence in 
some situations, the cavity field will reach a steady state with 
higher temperature than that for the incoherent case. 

Though we calculate the steady-state photon number in the 
cavity of the micromaser, we still emphasize that the motiva- 
tion of this paper is not to simply study the statistical prop- 
erties of the cavity field, but to study the quantum thermal- 
ization of a quantum system randomly coupled to a series 
of single-particle reservoirs in a time domain. Therefore our 
present work is different from other previous papers on quan- 
tum statistical properties of a micromaser (e.g., Refs.ll71-tl8ll). 



Here we employ the micromaser model only for convenience. 
The micromaser involves the process of a quantum system 
(the single-mode cavity field) randomly coupled with a series 
of single-particle reservoirs (these injected atoms). In other 
words, the micromaser model is the platform to show our idea 
of thermalization. More importantly, we focus on the temper- 
ature of the cavity field at a steady state. If the temperature of 
the cavity field at a steady state is equal to that of the injected 
atoms, we consider that the cavity field of the micromaser has 
been thermalized by these injected atoms. For the case of ran- 
dom injections, the steady state of the cavity field can be nat- 
urally identified as a thermal state with the same temperature 
as that of the injected atoms. 

In addition, from the viewpoint of experimental implemen- 
tation, this paper also provides a possibility for the examina- 
tion of thermodynamics with a cavity QED system. As we 
know, the cavity QED system has became a mature candidate 
for the implementation of experiments in quantum physics 
and quantum information processing Bill . Therefore the 
present work can also be considered an example for the exper- 
imental examination of thermodynamics with a cavity QED 
system. 

The paper is organized as follows. In Sec. [II] we present our 
thermalization model of a single-mode cavity field interacting 
with a series of atoms injected randomly. A quantum master 
equation is derived to describe the dynamics of the single- 
mode cavity field. In Sec. [TTT1 we show that the present quan- 
tum thermalization model can give a unified description of 
cooling, masering, and thermalization processes. In Sec. II VI 
we study the quantum thermalization when the initial state of 
the injected two-level systems possesses some quantum co- 
herence. In Sec.[V] we propose an experimental implemen- 
tation of our quantum thermalized model with superconduct- 
ing circuit-QED. We also show that the dynamics of the cav- 
ity field in the micromaser is equivalent to the dynamics of 
the transmission line resonator in the circuit QED. Finally, we 
conclude this paper with some discussions in Sec. [VI] 
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FIG. 1: (Color online), (a) Schematic diagram of our single particle 
thermalization model that a series of prepared TLSs randomly pass 
through a single-mode cavity one by one, in equilibrium it is equiv- 
alent to the conventional reservoir model, (b) where many identical 
atoms with spatially random distribution thermalize the single-mode 
cavity field. 



tially random distribution thermalize the single-mode cavity 
field. 

A single TLS interacting with the single-mode cavity field 
is described by the Jaynes-Cummings (JC) Hamiltonian 



H = "^~o"z + wo a + g(a& + + <r_a'), 



(1) 



where a and a' are, respectively, the annihilation and creation 
operators of the single-mode cavity field, they satisfy the usual 
bosonic commutation relation [a, a ' ] = 1. Hereafter we set 
ft = 1 . The operators of the TLS are defined as 



& + = frl = \e)(g\, & z = \e){e\ - \g){g\. 



(2) 



The parameter g is the coupling strength of the cavity field 
with a TLS. 

In the rotating picture with respect to 



Hq = -cr z + a>a a, 



(3) 



the Hamiltonian becomes 



II. CAVITY QED MODEL FOR THERMALIZATION WITH 
SINGLE-PARTICLE RESERVOIR 

The cavity QED model (as illustrated in Fig. [If a)) for ther- 
malization contains a single-mode cavity field of frequency u> 
and a series of injected two-level systems (TLSs) with excited 
state \e), ground state \g), and energy separation a>o- These 
TLSs pass through the single-mode cavity one by one ran- 
domly. Here, the single-mode cavity field is considered as the 
system to be thermalized, while the TLSs are considered as 
the single particle reservoir. The injections of the TLSs into 
the cavity are random and there is a limit of one TLS in the 
cavity each time. According to the viewpoint in statistical me- 
chanics, the average over an ensemble is equivalent to the time 
average in some sense. It is expected that the single-model 
cavity field will approach a steady state equilibrium with a 
temperature as that of the injected atoms, since this system is 
equivalent to the conventional thermalization model, as shown 
in Fig.QJb), where many identical atoms (reservoir) with spa- 



where 



Vi = 2°"* + g( ad ~+ + &- a )> 



6 = o)q — a> 



(4) 



(5) 



is the detuning of the cavity frequency a> with the energy sep- 
aration <dq of the TLS. In the resonant case, namely 5 = 0, the 
unitary evolution operator governed by the Hamiltonian (|4]i of 
the cavity QED reads ^ 



U(t) = exp(-/V / r) 



COS (gT Vflfl"'") - 



a 



-ia' 



cos (gT Vfltfl) 



(6) 



which is written in the Hilbert subspace of the TLS with the 
basis states 



\e) 



\g) 



(7) 
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To thermalize the cavity field, a series of TLSs are ran- 
domly injected into the cavity for a fixed time interval t. All 
of the TLSs are initially prepared in the density matrix 

Ptls = Pe\e)(e\ + p g \g)(g\ + A\e)(g\ + A*\g)(e\, (8) 

where A is the parameter describing the coherence of the 
TLSs. The state preparation of the TLSs can be realized by 
using a pumping field to excite the TLSs. We assume that the 
j'th TLS is injected into the cavity at time tj. After an interac- 



In fact, in addition to the action of the injected TLSs, 
the cavity inevitably couples with an external environment 
through the cavity wall. Within the quantum noise theory, we 
model the external environment of the cavity as a heat bath. 
When the coupling of the cavity field with the heat bath is 
weak, the decay of the cavity field can be described by lllOll 

_£p = —K(n,h + l)(2apa^ — a' ap - pa' a) 

+~Kn t h(2a i 'pa — acfip — pact*), (11) 

where k is the decay rate of the cavity. The thermal average 
photon number is 



with Pb = 1/(1(bT) being the inverse temperature of the heat 
bath. Hereafter we denote pit) as p to be concise. 

Since the TLSs are injected at random, we can introduce a 
rate r of a Poisson process to depict the arrival of the TLSs. In 
a time interval of (t, t + 5t), the probability of a TLS arrival is 
rdt. Hence the density matrix of the cavity field at time t + 5t 
can be written as lf32tl 

p{t + 6t) = (l- rSt) [pit) + £p(t)6t] + r5tMiT)pit). (13) 

Here the first term on the right-hand side of Eq. ( fT3l l describes 
the density matrix of the cavity field at time t + 5t when a TLS 
does not pass through the cavity, with the probability 1 - rdt. 
In this case, the cavity field evolves under the action of £.. 



tion of time t, the state of the cavity field becomes 

f>itj + T) = Tr TLS [UiT)p(tj)®pTLsUHT)] 

= M(r)p(tj), (9) 

where TriLS means tracing over the degree of freedom of the 
TLS. The superoperator Ai(j) introduced in Eq. (0 can be 
expressed as follows: 



(10) 



Additionally, the last term on the right-hand side of Eq. ( fT3l ) 
describes the density matrix of the cavity field at time t + 6t for 
the case of a TLS passing through the cavity, with the proba- 
bility rdt. Notice that here we approximately neglect the ac- 
tion of the heat bath on the cavity field during the process of 
the TLS passing through the cavity, since the time spent by 
each TLS in the cavity is assumed to be much shorter than the 
mean time between two injections of the TLSs. 

Taking the limit of 5t — > 0, we can obtain the following 
quantum master equation ll9T- [l6ll 

p = r (M(T) - l)p + £p (14) 

to describe the evolution of the single-mode cavity field. 



III. UNIFICATION OF COOLING, MASERING AND 
THERMALIZATION 

The evolution of the quantum state of the cavity field is 
governed by the master equation ( fT4l >. which depends on the 
initial state of the injected TLSs. Firstly, we consider the case 
where no coherence exists in the initial state of the TLSs, i.e., 
A = in Eq. ([8]). In the Fock state representation, the evolution 
equation for the diagonal elements P n = (n\p\n) of the density 



sin (gr Vafl 1 ') sin( gT' 

M(T)p(tj) = p e cos (gT Vfla 1 " jp(tj) cos (gT Vaa 1 " j + p e a^ — p(tj)- 



-a 

<ldd 1 Vflfl 1 
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+p g ■ — cip(tj)a) h pgCosygT ya^a) p(t j) cos (gT Va+a 1 
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+iA cos (gT Vflfl"'") p(tj)a' — iAa — p(tj) cos (gT Vfl f a) 

Vaa'" VaS 1 " 

sin(gTViat) sm(gT^/a¥) 
+z'/T cos (gT va'a)p(tj) a - iA ap(tj) cos (gT Vaa' J , 
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matrix p in the master equation ( TBI becomes 

P„ = rp e [cos 2 (gT Vn + l) P„ + sin 2 (gT Vn) P n -i] 
+rp g [cos 2 (gT Vn) P n + sin 2 (gT Vn + l) P„+i] 



-rP„ + [2«P„_! - 2(n + 1)P„] 



/f(n fA + 1) 



[2(n+ 1)P„ +1 -2nP„]. 



(15) 



Using the relation p e + p g — 1 and after some simple collec- 
tion, the above equation (fT3T > becomes 

P„ = -r sin 2 (grVn + l) (/? e P n - p g P n+ i) 
-k(h + 1) [h th P„ - (h th + l)P„+i] 
+r sin 2 (gT V«) (pA-l - P g P n ) 
+m [h, h P„-i - (h th + 1)-P«] ■ (16) 

The steady state solution P n — leads to the detailed balance 
condition and the relation 



r sin 2 (gT Vn) {p e P n -\ ~ PgPn) 
+kh [h 1h P„-\ - (n, h + l)P„] = 0. 



(17) 



Then the ratio R„ = P n /P n -\ between two neighboring photon 
number populations is obtained as 



R„ 



rp e sin 



[gT Vn) 



• Kn th n 



rp g sin 2 Vn) + K(n,h + l)n 



(18) 



We can understand such thermalization to the steady state with 
the definite population ratio ( fT8l ) as a temperature information 
transfer process from the TLSs to the cavity field, namely, the 
curve of -(lnR„)/(L> can explicitly reflect the information of 
the temperature of the TLSs. 

Such a temperature information transfer process can result 
in various coherent manipulations for quantum state engineer- 
ing. An example is the cooling of the cavity field as a gener- 
alized thermalization for all injected TLSs initially prepared 
in the ground state, i.e., p e — and p g — 1. In this case, the 
TLSs on the ground state will take away the energy of the cav- 
ity field and then cool it to reach a lower temperature defined 
by the decreased photon population 



h th l 



J (n t h + 1)1 + sin 2 [gr VI) r/K 



(19) 



where Po is determined by the normalization condition 
Tj7=qP" ~ 1- This generalized thermalization mechanism 
was even used to cool the nanomechanical resonator by the 
pulse-driven charge qubit [22]. Another example with p e = 1 
and Dg - shows the maser processes of the cavity field (0- 
H(||23|], which is represented by the amplified photon popula- 
tion 



™ sin 2 (gT V7) r/K + n th l 
-PoU- 



i=i 



(nth + 1)/ 



(20) 
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FIG. 2: (Color online). The ratio R„ versus the photon number n is 
plotted for gr = 0.05, 0.2, and 0.3. The injected TLSs are prepared 
in thermal equilibrium of inverse temperature f} - 2.898 (T = 200 



mK). Other parameters are set as w 
and n th = 0.008 (T b = 100 mK). 



1, k/oj = 10~\ r = 2 x 10- 



where Po is determined by the normalization condition 

In the absence of the cavity field dissipation, i.e., k — 0, 
Eq. (fL8b becomes R n = p e j p g , which is irrespective of both 
the index n and the average injection rate r. For this case, the 
temperature information of the TLS is perfectly transferred 
to the cavity field. For example, when the TLS is initially 
prepared in the thermal equilibrium with temperature T, that 



Pe(T) = 



exp(-y8tLt/2) 
2cosh(j6w/2)' 



P g (T) 



exp(J3cj/2) 
2 cosh(y0<x>/2) ' 



(21) 



where T = l/fe/J), then the population ratio R n = exp(-j6w) 
of the cavity field is independent of the index n, thus a thermal 
equilibrium has the same temperature T as that of the TLS. 

We give a physical explanation about the steady state of the 
cavity field in the absence of the cavity decay. When k — 0, the 
evolution of the cavity is governed by the following quantum 
master equation 

p = r (M(T) - \)p. (22) 
In the short r case, we can make the short time approximation, 

cos (gT VaS f ) S3 1 - (gT) 2 aa ] /2, (23a) 
cos (gT Vfl f a) 1 - (gT) 2 a [ a/2, (23b) 
sin (gT Vaa T ) sa gT Vao" 1 ". (23c) 
Up to the second order of r, the master equation (1221) becomes 
a Pg, 



2 

ap 
'r — 
2 



(2apa - pa' a — ci [ ap) 
(2a x pa—paa > — ao'p), 



(24) 



where a - r(gT) 2 . Now, the injected TLSs are prepared in a 
statistical mixture of the excited and ground states. From the 
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above equation (|24| |. we can see that the TLSs prepared in an 
excited state excite the cavity at an effective rate ap e , while the 
TLSs prepared in a ground state take away the energy excita- 
tion of the cavity field at an effective rate ap g ll33ll . By com- 
paring Eq. d24l > with Eq. (TT[ , we can see that the long time 
accumulation of the actions of the injected TLSs is equivalent 
to an effective heat bath with the inverse temperature 



co p g 



(25) 



Therefore, the cavity field can reach a steady state even in the 
absence of the cavity decay through the walls. It is of interest 
that the inverse temperature /?« of the effective heat bath of 
the cavity can be controlled by changing the populations p e 
and p g of the injected TLSs. 

In the presence of the cavity field dissipation, i.e., k + 0, 
generally, it is impossible to define a temperature for the cav- 
ity field in the steady state, since in this case the ratio R„ given 
by Eq. ( TT8b depends on n. In Fig. [2] we plot R„ versus photon 
number n for different gT. Clearly, for small gT, R„ shows the 
independence of the photon number n. Therefore, it is possi- 
ble to define an effective temperature for the cavity field when 
gT is small. 

In the short interaction time t limit, i.e., gT V» ^ 1 for 
all experimental accessible photon numbers n, we make an 
approximation sin 2 [gT -\/«) ~ (gT) 2 n, which results in an 
n-independent population ratio 



R = 



ap e + Kti th 
ap g + K(n, h + 1) 



(26) 



Thus for the TLS injection in thermal equilibrium, we can 
define an effective inverse temperature for the cavity field 



jSeff = 

which satisfies the relation 



1 



In R, 



mm{J3 b ,/3} < /3 eff < max{J3 b ,/3}. 



(27) 



(28) 



It means that the cavity field will approach a thermal equi- 
librium with an intermediate inverse temperature y8 e ff between 
those for the TLSs and the heat bath. Additionally, for the case 
of fib = P, the cavity field will approach a thermal equilibrium 
of yS e ff = /3. This result is reasonable from the viewpoint of 
quantum noise. A system coupled with two heat baths with 
different temperatures will reach an equilibrium with interme- 
diate temperatures between those of the two heat baths [34]. 



IV. QUANTUM COHERENCE ASSISTED 
THERMALIZATION 



In the above section, we study the generalized thermaliza- 
tion for the incoherent case, in which the injected TLSs do not 
possess quantum coherence. In this section, we study the co- 
herent case, i.e., A + 0. In this case, up to the second order of 
t the master equation ( fT4l can be reduced to 



in the short t limit, where the effective Hamiltonian reads 



eff 



gcV + fa, 



(30) 



with £ = rgTA. The superoperator J' is defined as 

STp = — yi(2a' pa - aa' p - paa') 

+ —y2(2apa'—a'ap — pa^a), (31) 

where we introduced two transition rates: the decay rate j\ 
and the excitation rate 72, 



J\ = ap e + KHth, 

y 2 = ap g +K(h th + 1). 



(32a) 
(32b) 



During the derivation of the master equation d29l ), we have 
used the approximation given in Eq. ( l23l . 

The above effective Hamiltonian H e g describes the role of 
the quantum coherence of the injected TLSs: the off-diagonal 
terms in the initial state offer nonvanishing atomic transition, 
which is added as a driving source of the cavity field. The gen- 
eralized master equation d29l ) describes a driven cavity field in 
contact with an effective bath characterized by two rates. This 
effective bath consists of a TLS reservoir and a heat bath. It 
is worth pointing out that the properties of the TLS reservoir 
can be manipulated through changing the initial populations 
p g and p e . 

From the master equation d29l ), we can obtain the follow- 
ing equation of motion of the average value of the creation, 
annihilation, and photon number operators, 



dt 
d 



(a(t)) 



(72 - 7i) (0> - i£ 



dt 
d 

dt 



<a T (0) = ~ (72-70 <a T (*)> + if, 



(33a) 
(33b) 



(h(t)) = - (7 2 - y0 (ft (?)> - i?<a f (?)> + if {a (?)> + 7i ■ 

(33c) 



The steady state solutions of the above equation are 

2^ 



(a) ss = (iVy ss = -■ 



(n) ss = 



72 -7i 



(72 - 7i ) 72 ~ 7i 



(34a) 
(34b) 



p~i[p, H eff ] + Jp 



(29) 



The above results show that the quantum coherence can in- 
crease the steady state average photon number (with the first 
term in Eq. d34bb ) in the cavity field. When no TLS is in- 
jected, i.e., r — 0, the average photon number at steady 
state (n) ss = n t h- On the other hand, for the case of ther- 
mal TLSs injection of inverse temperature ft and a perfect 
cavity, i.e., k — 0, the steady state average photon number 
(ft) ss = l/[ ex P(jSw) - 1], which implies that the cavity ap- 
proaches an equilibrium with the same temperature of the 
TLSs. 

The above argument based on the short time approxima- 
tion is only a heuristic analysis, thus we need to numerically 
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FIG. 3: (Color online). The -(lnP„)/a> versus the photon number 
n is plotted for both the coherent (A = 0.08) and incoherent (A = 0) 
cases. Here, gr = 0.05 and other parameters are set as those in Fig. [2] 
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FIG. 4: (Color online), (a) Schematic diagram of the circuit QED 
system of a transmission line resonator (TLR) coupled with a charge 
qubit, which is controlled through an external magnetic flux O and a 
gate voltage V. (b) The energy levels of the charge qubit versus the 
magnetic flux <t>, the working status of the qubit is controlled by the 
magnetic flux "pulse series". By tuning the magnetic flux O, we can 
switch on or off the coupling between the TLR and the charge qubit. 



solve the master equation < fl~4T > directly. In Fig. [3] we plot the 
-(lnP„)/u> versus the photon number n for the coherence and 
incoherence cases. It can be seen from Fig. [3] that the coher- 
ence in the initial state of the TLSs can increase the steady 
state temperature of the cavity. However, for the coherence 
case, the dependence of -(In P„)lu) on the photon number n is 
approximately linear therefore it is an approximation to define 
an effective temperature for the cavity field. 



V. EXPERIMENTAL IMPLEMENTATION WITH CIRCUIT 
QED SYSTEM 

In this section, we present an experimental implementation 
of our quantum thermalization based on superconducting cir- 
cuit QED system 1I35I l36Tl . As shown in Fig.Ua), a supercon- 
ducting transmission line resonator (TLR) couples to a super- 
conducting charge qubit. After the quantization of the elec- 
tromagnetic field in the TLR, the Hamiltonian describing a 
single-model field in the TLR reads lf35ll 



H 



TLR 



(35) 



where a> is the resonant frequency of this mode. Here we 
only choose the single mode which is (near) resonant with 
the lowest two levels of the superconducting Cooper-pair box, 
i.e., the charge qubit. The Hamiltonian of the Cooper-pair box 
(CPB) superconducting circuit is IB7I1 

#cpb = 4E c Y(n-n g ) 2 \n){n\-EjCOs(^-\ 

x2(|« + l><«| + |»X» + l|), (36) 

neZ 

where E c = e 2 /(2C£) is the Coulomb energy, with Cy, = 
2Cj + C g being the total capacitance connected with the su- 
perconducting island. Ej is the Josephson coupling energy 
of a single Josephson junction. An external magnetic flux 
O through the loop can tune the effective Josephson cou- 
pling energy of the CPB. The symbol Oq is the magnetic flux 



quanta. In addition, we introduce the gate Cooper-pair num- 
ber n g = C g V/(2e) with the gate capacitance C g and the gate 
voltage V. The state \n) (n e Z, where Z denotes the integer 
set) stands for n extra Cooper pairs on the island. 

In the present scheme, since the existence of the TLR, the 
gate voltage contains two parts: a dc part V g c and a quantum 
part Vq generated by the TLR. Then 



2e 



— n„ -\ , 

* 2e 



(37) 



where n d g c = C g Vf'/(2e). The voltage generated by the TLR 
can be written as 



(38) 



Lc 



where L is the length of the TLR and c is the capacitance per 
unit length of the TLR. 

With the substitution of Eqs. (ITTl i and ( 1381 into the Hamil- 
tonian (f36t and restriction into the subspace with basis states 
|0) and |1), we obtain 

Hcpb = -2E C {\ -2nf)f z -£ y cos^Jf, 



Cy V Lc 



(fl + fl T )(l -2nf -r z ), (39) 



where we have introduced the Pauli operators 

^ = 11X01 + 10X11, f z = |0><0| -|1><1|, (40) 

and discarded some constant terms. When the charge qubit is 
working at the optimal point n d g =1/2, the total Hamiltonian 
of the system becomes 

t /7rO\ C g I cj + 

H = u>a'a - Ejcos — — t c + e— - ,/ — (a + a')T z . 

\ 0>o / C 2 V Lc 

(41) 
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By making a rotation 

f x — > — fr z 



the Hamiltonian becomes [351 

/i = wo 'a + — cr : + g(a + a')cr x . 
where the energy separation of the charge qubit is 



n<t> 

too = 2Ej cos — 
and the coupling strength is 



,5k 
c 2 



(42) 



(43) 



(44) 



(45) 



The Hamiltonian d43T > reduces to the usual JC Hamiltonian 
given in Eq. (Q]) by making the rotation wave approximation. 

From Eqs. (144-b and ((45), we can see that the energy separa- 
tion w of the charge qubit is tunable by controlling the bias- 
ing magnetic flux <£, and the coupling strength g is fixed once 
the superconducting circuit is fabricated, therefore the effec- 
tive method to switch on and off the coupling of the resonator 
with the charge qubit is to tune <t> such that the qubit cou- 
ples with the resonator in resonant and very largely detuned, 
respectively. We schematically plot the qubit energy level ver- 
sus the magnetic flux <t> in Fig.HJb). The working points "on" 
and "off" correspond respectively to the resonant coupling and 
decoupling of the qubit with the cavity. The magnetic flux 
"pulse serial" controls the interaction of the qubit with the res- 
onator. Note that similar methods have been proposed to gen- 
erate photon Fock states ID8I1 and have recently been realized 
based on a circuit QED system consisting of a transmission 
line resonator coupled with a phase qubit ll39l l40Tl . 

However, it should be emphasized that, strictly speaking, 
the present circuit QED system is different from the mi- 
cromaser system since a micromaser has many independent 
atoms, while the TLR has strictly one TLS. Therefore we 
need to know the conditions under which the dynamics of the 
single-mode field in the TLR is equivalent to that of the cav- 
ity field in the micromaser. Without loss of generality, in the 
following we study the dynamics of the circuit QED during a 
single cycle. We assume that the coupling of the charge qubit 
with the TLR is switched on at time f,, and after an interaction 
of time t, this coupling is switched off. We denote the density 
matrix of the circuit QED at time f, + t as 



Pcir-QED = Pmnrs I'«>TLR<«ItLR 1 

m,n-Q r,s={e,g) 



\r)M 



(46) 



where states |;«(«))tlr and \e(g)) q denote the states of the TLR 
and the charge qubit, respectively. Generally, this density ma- 
trix given in Eq. ( 1471 ) is an entangled state due to the coupling 
between the charge qubit and the TLR. During the time in- 
terval from f, + r to ti+u the time of the (i + l)th turning on 
the coupling, the density matrix (f4Tb evolves under the local 
actions of the environments of the charge qubit and the TLR. 



In a micromaser, correspondingly, we only focus on the 
quantum state of the cavity by tracing over the atom. There- 
fore, if the total density matrix of the cavity and the ith in- 
jected atom is 



Pcav-QED 



z z 

m,n=0 r,s={e,g 



:l'M>cav<«lcav ® |r) aj (4j, (47) 



at time ?,■ + r, where |m(n)) cav and \e{g)) a . denote the states of 
the cavity and the zth injected atom in the cavity QED, respec- 
tively, then the reduced density matrix of the cavity at time 
ti+i should be 



Tr fl; [/)cav-QED] 



^ ^ ^ ^ Pmnrrl'^)cav('^lcav 

m,«-0 r={e,g\ 



(48) 



taking the trace over the z'th injected atom. Notice that where 
we have neglected the action from the environment of the cav- 
ity during the interval from f,+T to f,- + i . In addition, the (i+l)th 
atom is prepared in its initial state at time fj+i. Therefore, to 
simulate the micromaser with the circuit QED system, it is re- 
quired that the qubit should be disentanglement from the res- 
onator at time f,-+i to avoid the correlation between the qubit 
and the resonator at the beginning of the (i + l)th coupling. 
By comparing Eq. d47T i with Eq. d48l i. we can see that the dis- 
entanglement condition is that, during the time interval from 
t i + t to a time f,- +t+ r r before f, + 1 , the qubit should be relaxed 
to its ground state as follows: 



\e) q {e\ q -» \g) q (g\ q , 
\e) q (g\ q -» 0, 
\g) q {e\ q -» 0, 

\g)q(g\q -* \g)q(g\q- 

Under this process the density matrix d47l > becomes 



Pcir-QED = PTLR ® \g)q(g\q, 

where we denote 

oo 

PTLR = ^ ^ Pmnrr\m)TLR{n\jLR- 

m,n=0 r— {e,g} 



(49) 



(50) 



(51) 



Clearly, the density matrices ptlr of the TLR in the circuit 
QED and p cav of the cavity in the cavity QED have the same 
form. 

Before the beginning time of the (i + l)th coupling, we 
need to prepare the qubit in its initial state, we denote the state- 
preparation time is t p . At time f, + i, we switch on the coupling 
between the qubit and the TLR, repeating the process as de- 
scribed before, we can simulate the micromaser with the cir- 
cuit QED system. 

Now, there are several time scales in the circuit QED sys- 
tem: the interaction time r, the relaxation time r r , and the 
state-preparation time x p . From the above discussions we can 
see that the requirement of these time scales is 



T + T r + T p < t i+ l - ti 



(52) 
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for all i. In addition, since we neglect the relaxation of the 
qubit during the interaction time r, then the relaxation time 
T r » t. The preparation of the qubit' s state can be realized 
by using a classical field. For example, a 2n pulse can trans- 
fer the qubit from its ground state \g) to excited state \e). The 
preparation time t p can be much smaller than other time scales 
through choosing a sufficiently strong field. Namely, we ap- 
proximately have the relation t + t,- + t p ~ t + r r . 

In the following, we give a simple estimation of the above 
time scales under the current experimental conditions. Ac- 
cording to recent circuit QED experiments B36I1 . we take the 
following parameters, the resonator frequency u> - 2n x 10 
GHz, the coupling strength g = 2n x 50 MHz, the cavity de- 
cay rate k — 2n x 1 MHz. The interaction time r ~ 10~ 9 s 
for gT ~ 0.05. The rate to switch on the coupling r ~ 2n x 2 
MHz. Namely, the interaction between the resonator and the 
qubit takes place every t* 1/r ~ 10~ 7 s. In other words, the 
average time f, + i - f, « 10~ 7 s, where the overline represents 
average value. Then the qubit relaxation time r r should be of 
the order of 10 ns to satisfy the conditions t + t> < f J+ i - f, 
and T r » t. Therefore, the requirements for the design of 
the random pulse are min{f !+1 - f,} > t + r r ~ 10~ 8 s and 
ti+i - t, ~ 10~ 7 s. We can decrease the average rate r to en- 
large the variable space of the time scale r r . As an example, 
we take the temperatures of the TLS and the field heat bath as 
T = 200 mK (J3 = 2.898) and T b = 100 mK (fi b = 4.797), re- 
spectively. We can calculate the temperature of the cavity field 
at thermal equilibrium as T e g = 100.54 mK (J3 e s = 4.773) for 
the case gr = 0.05. 

VL CONCLUSIONS AND DISCUSSIONS 

In conclusion, we have proposed a kind of generalized ther- 
malization with a single-particle reservoir. This generalized 
thermalization uniquely describes the cooling, masering, and 
thermalization processes. We have shown our generalized 
thermalization based on a micromaser-like system, in which 
a series of well-prepared TLSs are injected randomly through 
the cavity. In the absence of the cavity decay, the cavity can 
reach an equilibrium with the same temperature as that of the 
TLSs. When the cavity is coupled with a heat bath, at a steady 
state the cavity can reach a thermal equilibrium with an inter- 
mediate temperature between those of the heat bath and the 



TLSs. We have also studied the effect of quantum coherence 
on the thermalization. It was found that the quantum coher- 
ence can increase the temperature of the system to be ther- 
malized at a steady state. We have suggested an experimental 
implementation of our generalized thermalization with the su- 
perconducting circuit QED system. 

We point out that the present investigations have some po- 
tential values in solid thermodynamical applications. For ex- 
ample, in circuit QED we can manipulate the steady state of 
the TLR by preparing the quantum state of the qubit. When 
the qubit is prepared in its ground state, the steady state of 
the TLR will be in its ground state. A similar idea has been 
used to cool a TLR 114 111 . In addition, by preparing the qubit 
in its excited state, the TLR can be used to realize a solid laser 
on-chip 10]. 

We give some discussions concerning the ergodicity during 
the quantum thermalization process. As mentioned in the In- 
troduction, one of the motivations of our present investigation 
was to show the quantum thermalization of a system randomly 
coupled with a series of single-particle reservoirs in a time 
domain. Essentially, the physical principle at the background 
is the ergodicity, which involves the equivalence between the 
time average and ensemble average j43j. For the present case, 
it has been shown that the form of the quantum master equa- 
tion (2M obtained in the short t limit can guarantee the er- 
godicity I44I l45ll . Therefore, it is a correct result that these 
injected atoms can thermalize the single-mode cavity field in 
the micromaser model. 

Finally, we emphasize that in this work the injected atoms 
are prepared in thermal equilibrium with positive temperature. 
As is known, the temperature of two-level atoms can be nega- 
tive 14611 . However, for a harmonic oscillator, it is impossible 
to define a negative temperature since the energy of the har- 
monic oscillator is finite. Therefore, the discussions in this pa- 
per are restricted within the case that the population p e of the 
excited state is smaller than the population p g of the ground 
state. 
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